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i ' Abstract 
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We outline the existing descriptions of the charm component of the deep inelastic 
proton structure function F2. We discuss recent approaches to include charm mass effects 
^ ; in the parton evolution equations and the coefficient functions. 
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Both HI and ZEUS have used c — > D* to measure charm production in the HERA domain 
They find F£/F 2 ~ 25%, as compared to the EMC measurement of F^/Fc ~ 1% at 
lower fixed-target energies. Moreover, the precision of the measurements of F£(x, Q 2 ) will 
improve, particularly when vertex detectors become available. Such measurements should reveal 
important information on the gluon, which enters at LO through the photon-gluon fusion (PGF) 
process 7$ — > cc. Clearly it is important to see how best to include the charm mass effects in 
the analysis. At NLO we have 



FZ(x,Q 2 ) = I' dz- \lc q=c {z,Q\^)c(-,A + \c g {z,Q\f)g(-,^ 
j x z i_y \z j y \z 



(1) 



where the coefficient functions C c = CfP' + asC^ and C„ 



(a) Massless charm evolution 



The simplest approach is to assume that the charm distribution c(x,fi 2 ) = for /i 2 < 
where fi c ~ m c , and then to evolve assuming c is a massless parton both in the splitting and 
coefficient functions. This is the approach used in the MRS and all but the latest (CTEQ4HQ 
p|) of the CTEQ global analyses. MRS take /x 2 = 2.7GeV 2 so as to obtain a satisfactory 
description of F 2 C of EMC ||. On the other hand CTEQ set /i c = m c ; this is a consequence of 
their choice to use the ACOT scheme 0, ||, explained below, to define the parton densities. In 
both cases, the number of quarks in the hard scattering coefficients is the number of "active" 
quarks. 

Corresponding conditions are applied for the b quark. This is a minor issue in the global 
fitting because of the small contribution of processes involving b quarks. 

Although phenomenologically successful, the massless model clearly is inadequate in the 
charm threshold region. For instance an on-shell cc pair can be created by photon-gluon fusion 
(PGF) provided 

W 2 = Q 2 (l-x)/x>4m 2 c 

where W is the 7*g — > cc centre-of-mass energy. Thus at small x, cc production is not forbid- 
den even for Q 2 < [i 2 where the massless approach gives zero. The physical threshold for cc 
production, W 2 = ^m 2 c1 is not the threshold that is provided by the massless model. Only to 
the extent that charm production is small in the threshold region does the massless model give 
a useful approximation to F 2 . 



(b) Photon-gluon fusion 



In this approach charm is treated as a heavy quark and not a parton. That is we put c = 
for all Q 2 in (P and use the known m c ^ gluon coefficient function Cj GF for C g . This is called 
a fixed flavour number scheme (FFNS) with n./ = 3. In contrast to the situation for massless 
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quarks, there is no collinear divergence in 7*^ — > cc since the integral over the cc transverse 
momentum is regulated by m c . However, this in turn means that at large Q there are large 
logarithmic corrections in higher orders: F 2 ~ [as ln(Q 2 /m 2 )] n at n th order. In terms of the 
FFNS all these contributions should be absorbed in the coefficient function. However, these 
are just the large logarithms that are resummecQby DGLAP evolution in the purely massless 
approach. This implies that at large Q 2 we should treat charm as a parton. 

(c) Variable flavour number schemes: ACOT and MRRS 

The aim is to obtain a universal charm distribution c(x, Q 2 ) all the way down to the reso- 
lution threshold Q 2 ~ m 2 in order to make predictions for other processes. To do this we must 
include the nonzero mass of the charm quark in the calculation in a consistent way. 

The ACOT [|], |JJ method is to exploit the techniques of Collins, Wilczek and Zee ||. Below 
fi = m c , the parton densities are those of the FFNS with n/ = 3, and above /1 = m c , they are 
pure MS distributions.^] The evolution coefficients are independent of mass in the MS scheme, 
so they are the same as in a purely massless calculation, with the number of active flavors 
changing from 3 to 4 at jj, = m c . At leading order explicit calculation |4| shows that the parton 
densities are continuous at \i — m c . The order a 2 s corrections to the matching conditions have 
been calculated J7| but not yet implemented. 

In the ACOT scheme, the remaining dependence on the nonzero charm-quark mass appears 
in the coefficient functions. In the threshold region, the charm density c(x, \i — Q) by itself 
does not accurately provide the charm contribution to the structure function F£ . The necessary 
correction, at the first non-trivial order is provided by the gluon coefficient in Eq. |l[ Consistency 
with calculations in the FFNS is obtained at the appropriate level of accuracy. 

Buza et al. |7j have evaluated the coefficient functions for Q 2 ^> m 2 in the VFNS and the 
FFNS to NLO. Their definition of the parton densities is exactly that of ACOT ||. They find 
that for Q 2 > 20GeV 2 the VFN and FFN schemes agree very closely. They also work out the 
matrix relation in the limit Q 2 3> m 2 which allows a matching of all the parton four flavour 
densities above the scale /1 2 = m 2 , to the three flavour densities below that point. A direct 
calculation of the matching conditions |H| from the partonic operator matrix elements confirms 
the matching conditions computed by Buza et al., for the subset of the cases that have so far 
been computed. This relation may be used to calculate other inclusive cross sections (e.g. large 
Et jets etc.) This particular calculation only applies for Q 2 3> m 2 , which does not solve the 

*The NLO corrections to the PGF structure function are known || and a leading twist analysis has been 
used to perform a resummation of the [as \n(Q 2 /m 2 )] n terms for Q 2 3> m 2 c . Unfortunately such an approach is 
not applicable for F2 in the threshold region Q 2 ~ m 2 . 

t Note that the definition includes the full unapproximated dependence of the parton densities on the charm 
mass, and that "MS" refers to the ultra-violet renormalization of the parton densities. Also, the position of the 
change of definition need not be at /j, — m c ; any other value in the neighborhood would be suitable, provided 
that corresponding changes in the matching conditions are made. 
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problem of what to do at Q 2 ~ m 2 . One of the authors (JCC) considers the principles of 



the problem solved; a paper is in preparation [O]. Work is in progress to obtain the ACOT 



coefficient functions at order a 2 s from the Buza et al. calculations; these coefficient functions 
will apply at all values of Q of order m c and larger. Another author (MGR) considers that there 
has to be another matching condition which provides parton densities which are continuous in 
the threshold region. 

In common with the coefficient functions used in all the global analyses, those of ACOT are 
for inclusive cross sections, F 2 for the case under discussion here. So a different, but related 
calculation must be made if, for example, one wishes to compute the distribution of charm 
in the final state. f\ The reason why it is necessary to treat each process separately can be 
illustrated by this example. The charm pt spectrum associated with the first cc pair produced 
by evolution from a gluon is proportional to 



dp T 



Pt + m 1 . 



2 , 2 2ml z(l - z) 

z z + (1 - z) + 



pip + 



(2) 



This same pt spectrum appears in the calculation of the ACOT matching conditions for the 
parton densities (integrated over p?). After application of the rules for computing an MS 
distribution, ACOT obtains the matching condition given previously. However only the first 
term in eq. (§) is operative in giving the (mass-independent) evolution kernel of the ACOT 
distributions. The second term is hidden in the matching conditions and coefficient functions. 
The full contribution to F 2 C at order a s also needs the gluonic coefficient function in eq. ([[]). 

The problem that worries MGR is that the second term in (|2]) is proportional to as{rri 2 ). 
Thus it cannot be absorbed in a fixed order coefficient function, which necessarily depends on 
o<-s{Q 2 )- On the other hand the matching takes place at one fixed point Q 2 = m 2 , where we 
deal with as(m 2 ). It is possible to choose such a matching which gives the correct description 
for Q 2 ^> m 2 but it looks impossible to correct the whole pt spectrum in the threshold region 
by just changing the parton densities at one point. 



A different approach to include m c ^ effects has been proposed by MRRS |L2] . Their aim 
is to formulate the evolution procedure with m c ^ so as to generate universal parton distribu- 
tions applicable to any inclusive or exclusive process with the (known) coefficient functions in 
the conventional MS renormalisation scheme. In order to do this they analyse the leading log 
contributions which come from the relevant (ladder type) Feynman diagrams. It turns out that 
a remarkable simplification occurs. Recall that the leading logs come from the configuration 
where the pt of the partons are strongly ordered along the ladder. At NLO accuracy it was 

* A closely related case is that of the Drell-Yan cross-section, da^y / dQ 2 dy integrated over transverse mo- 
mentum qT- The physical cross section is analytic at qp = 0, but the lowest-order, or parton-model, contribution 
is proportional to <5(<7t)- Higher order terms give a cross section da/d^q that, without resummation, diverges 
as qr — > and that has delta-function with divergent coefficients. Integration over qT cancels the divergences 
and gives a correct estimate of da/dQ 2 dy. 
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found to be sufficient to take into account the m c ^ effect for the charm parton with px ~ m c 
In fact all the m c ^ effects at NLO occur onlyfj] in P cg 



p 

1 eg 



z 2 + (l-z) 2 + 2 -^z(l-z) 



(3) 



with Q 2 = ml +p T . In this way the MRRS procedure automatically reproduces the correct px 
spectrum of charm quark. It is straightforward to generalise the MRRS to higher orders ||12|| . 



In a VFNS we have to take care to avoid double counting. The same Feynman graphs are 
generated by the evolution with zero-order coefficient function (corresponding to 7*c — > c with 
a spectator c, or vice versa with c <-> c) and the first order coefficient function (which describes 
photon-gluon fusion (g — > cc) ® (7*0 — > c)). Thus we must subtract the contribution generated 
by evolution from the PGF contribution. This is done by both ACOT and MRRS. 

The charm component F 2 C of F2 is totally determined by the gluon and light quark densities. 
The only parameter is the value of m c . The MRRS predictions are shown in Fig. 1. They are 
in good agreement with both HERA and fixed-target EMC data, which together cover a wide 
range of x and Q 2 . Unfortunately in order to match on to conventional MS coefficient functions 
the evolution scale is required to be Q 2 = ml + p\ leading to a rather low charm threshold at 
Q 2 = ml. The value of as is not sufficiently small in this region to neglect NNLO corrections. 
One of the main effects is to move the threshold to Q 2 > Am 2 c which, due to the Uncertainty 
Principle, is the virtuality of the photon required to resolve a charm quark within the g <-» cc 
fluctuations. In the NLO predictions in Fig. 1 the NNLO modification has been imposed by 
hand in a simplified way. This is the origin of the artificial structure at Q 2 = Ami. Now that 
a complete NLLO calculation is available it may be used to smooth out this behaviour. 

A corresponding plot from CTEQ, with the use of the ACOT scheme, is shown in Fig. 2. 
It compares the experimental data on the charm contribution to DIS with calculations with 
two sets of parton densities. One set, CTEQ4M, is obtained from a global fit which uses the 
massless scheme, and the other, CTEQ4HQ1, is obtained from a global fit [0 that used the 
ACOT scheme for DIS. 



§Apart from the corresponding adjustment in the 6(1 
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— z) term in P gg JO. 
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Fig. 1: The predictions of MRRS @ for F 2 C compared with the EMC @ and HERA (T| measurements. 
The dotted, continuous and dashed lines correspond to m c = 1.2, 1.35 and 1.5 GeV respectively. The 
starred data points in the HERA domain are obtained by interpolating ZEUS measurements, whereas 
the other HERA data correspond to the HI measurements. 
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Fig. 2: The CTEQ predictions [14] for F| when the ACOT scheme is used. The dashed line gives the 
prediction from the CTEQ4M distributions, while the continuous line shows the prediction from the 
CTEQ4HQ1 distributions, which were obtained from a global fit Q where the ACOT method was 
used. The data is from HI ill. 
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